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Abstract 


X dimensionless axial coordinate, x*/h* 


Sound propagation without flow in a rectan- 
gular duct with a converging-diverging area varia- 
tion is studied experimentally and theoretically. 
The area variation was of sufficient magnitude to 
produce large reflections and induce modal scatter- 
ing. The rnts (root-mean-squared) pressure and 
phase angle on both the flat and curved surface 
were measured and tabulated. The "steady" state 
finite elenient theory of Astley and Eversman and 
the transient finite difference theory of White are 
in good agreement with the data. It is concluded 
that numerical finite difference and finite element 
theories appear ideally suited for handling duct 
propagation problems wnich encounter large area 
variations. 
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List of Symbols 
speed of sound, m/s 
frequency, Hz 

duct height ahead of area change, m 
maximum height of area change, ni 
function of x and y. Ref. 12 
dimensionless reduced frequency, 

iu*n*/co 

length of area variation, m 


y 

Y' 

Ze 

«r 


a.a.Yl'l 
o,T J 


* 

Po 


dimensionless transverse coordinate, y*/h* 

height of area variation, Eq. (1) 

exit acoustic impedance, kg/m^sec 

sound power reflection coefficient, 

Eq. (10) 

function of x and y. Ref. 12 


specific heat ratio 

transverse coordinate in transformed 
plane, see Fig. 5 

axial coordinate in transformed plane, 
see Fig. 5 


reference density, kg/m^ 
ui dimensionless frequency h*f*/co 

ui* angular frequency, rad/ sec 

Superscript; 

* dimensionless quantity 

Introduction 


P time dependent dimensionless pressure, 

P*(x*,y*.t*)/pS cJ2 

p Fourier transform pressure, 

P*(x,y)/po C$2 

Po mean pressure, p$(x,y)/pg cg 2 

S standi DQ wave ratio 

t dimensionless time, t*f* 

Uq mean axial velocity, Uq/cS 

u acoustic axial velocity, u*/cq 

Vq mean transverse velocity, Vq/co 

V acoustic transverse velocity, v*/co 


“Steady" state and transient finite difference 
and finite element theories have been developed to 
study sound propagation in complex ducts with axial 
variations in cross-sectional area, wall liner 
impedance (absorbers) and with gradients in flow 
Mach number. Reference 1 contains a comprehensive 
deSL.T lotion cf the techniques, advantages, and 
limitations associated with the various numerical 
solutions of the sound propagation equations. In 
general, the literature is concerned primarily with 
theoretical solutions of the sound propagation 
equations. The present paper concentrates on the 
experimental verification of the numerical theories 
particularly for ducts with large area variations. 

Only a limited number of exf>erimental studies 
is available to verify the numerical as well as tne 
analytical sound propagation theories. For 
straight soft wall ducts, measured transmission 
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losses and axial. pressure variations are in reason^ 
able agreement with analyticai^*^ and finite 
element^ theories. For acoustic filU/s and 
mufflers with abrupt area changes, finite element 
predictions of transmission losses^t® are also in 
good agreement with theory. Kagawa, et a1.' 
recently developed a combination of the finite ele- 
ment (In duct) and analytical methods (Greens- 
Theorem In the far field) to analyze sound propaga- 
tion from a loud speaker. The measured far field 
sound pressure and directivity characteristics were 
predicted within a few decibels over a wide fre- 
quency range. In a much simpler approach, Ref. 8 
relates In duct finite element solutions to trends 
in the radiated far field data. 

Nayfeh, et al.^^10 experimentally studied 
sound propagation in an annular duct with a center- 
body with a slowly varying cross sectional area. 
This experimental test section was designed to corh 
ply with the assumption Inherent in their multiple 
scale perturbation analysis. The centerbody slope 
was limited to 0,2 to reduce reflection? and 
coupling among acoustic modes, both of '•'^hich are 
neglected In the analysis. The circumferential 
variations of pressure amplitude and phase at 
several axial positions were examined with bard 
walls and lined sections with and without mean 
flow, Reasonable agreement between the theoretical 
and experimental results were found for the con- 
verging portions of the duct. 

The present study involves a similar experi- 
ment with axial area variations buc ^ hard wall 
rectangular duct with no flow. In contrast to 
Refs. U and 10, however, large slopes In the duct 
wall will be introduced to enhance reflections and 
modal scattering to obtain a severe test of the 
theory. Also axial pressure measurements are made 
in both the converging and diverging portion of the 
duct* 

In the present paper, the experimental test 
section, apparatus and procedure will first be dis- 
cussed. The experimental data will be conveniently 
tabulated for later use In the literature. Then, 
one steady state finite element thepry^^ and one 
transient finite difference theory^^^ will be com- 
pared to the data and conclusions drawn. 

Experimental Apparatus and Procedure 

To verify the numerical theories for hard 
walls without mean flow, the general acoustic duct 
system descrioed in^'^ was modified into the sIiik 
pie no flow apparatus shown in Fig. i. Tne system 
shown in Fig. 1 was aoapted to simulate plane wave 
propagation in an infinite duct. The following 
paragraphs describe the acoustic system including 
components and general operation. 

Source 

Sound was generated by a 120 watt, 300 Hz to 
6 kHz driver. In the experiment reported herein, 
the sine wave generator provides a 1560 Hz 

signal to the amp1if'>er which in turn drives the 
speaker shown in fig. 1. The duct temperature was 
at 27“ centigrade. Since the first nonplane mode 
begins propagating at a frequency of 1701 Hz, the 


choice of 1560 Hz for the driving frequency guaran- 
tees that only plane waves will propagate in the 
straight portion of the duct far from the area 
variation. 

Test Section 

The sound travels through a 10 cm by 3.81 cm 
(4 in. by 1.5 in.) rectangular test section. This 
1,4 meter (56 in.) long section consists of Id flat 
detachable plates (8 on top and 8 on the bottom). 

The variable area test piece, to be described short- 
ly, was Inserted in place of the bottom piece at 
location 6, An exponential horn was attached to the 
upstrean, end of the test section. For these experi- 
ments, the horn was fitted with an acoustic foam 
wedge to approximate a ppCo termination, which 
simulates an Infinitely Tong duct. 

Variable Area Test Piece 

The variable area test piece shown in Fig. 2 
was constructed from wood. The surface profile is 
prescribed by a fourth degree polynominal of the 
form: 

Where the length L* is 7.62 cm (3 In.) with a 
height H* of 1.91 cm (0.76 in.) which Is just half 
the height of the duct. This curve has zero slope 
at x*/L* equals 0, .6 and 1, and has a maximum 
slope of 0.77 at x*/L* equals 0.21 and 0J9# The 
symbols used here are defined in the list of symbols. 

This test piece has provision for measuring the 
pressure along the curved surface at x*/L* equals 
0., 0.3, 0.5, 0.7, and 1,0. The microphone holes 
at x*/L* equals 0.3 and 0,7 correspond wUh tne 
position of the peak component of the transverse 
pressure profile, as displaye-i in Fig. 5 of 
Ref, 12. The holes were stopped with wooden plugs 
when the microphones were not in place. 

i^icrophone Installation 

T.wo 0.64 centimeter (1/4 in.) diameter corn* 
denser microphones were used to determine the acous- 
tic field. The microphone designated A was used to 
monitor tfia source strength of the signal while 
microphone B was used to determine the pressure and 
phase angle along both the variable area test piece 
and along the flat surface above the variable area 
test piece, Tne output of each microphone was 
analyzed to give the rms (root-mean-squared) pres- 
sure and phase angle. The flat microphone test 
holder plate is shown in Fig. 3. This holder Is 
mounted above the variable area test piece in loca- 
tion 6 of Fig, 1. This aluminum plate was fabri- 
cated with 0.64 cm (1/4 In.) diameter holes which 
allowed one microphone (b) to be traversed in the 
axial direction in increments of 0.953 cm 
(3/8 in.). The test piece was fabricated with an 
insulator in each hole, so that the microphones were 
electrically insulated from the test plate. During 
a measurement, each microphone’s membrane was posi- 
tioned flush with the surface of the test plate 
(protective shields were removed). 

The response of these microphones was flat to 
25 kilohertz. A standard acoustic source was used 
to check the system amplification. This source 
generated a 114 decibel signal In 5 octave steps 
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from O.IiSS to 2.0 kilohertz. The microphone was 
placed In the source and the resulting decibel level 
of the output was noted. In this way an absolute 
scale wjs established for the Sf»L. In this experi- 
ment, the SPL level was set at 100 dB at a position 
0.9&3 cm (3/B In.) In front of the start of the 
variable area cross section. The microphone used to 
monitor the source has been designated as m1cro< 
phone A. 

General Procedure 


The equipment shown In Fig. 1 was turned on 
approximately four hours before data were taken. 

The rms pressure and phase angle were then deter- 
mined by microphone B In sequence at each axial 
location along both the flat and curved surfaces. 
Three measurements were averaged to determine a data 
set. As a further check on the repeatability of 
vhe measurements, a second set of measurements was 
recorded (data set II). The data Is displayed In 
Tables I and II, 

finally, the data In Table I have been arbi- 
trarily normalized to an average value of 1.4. This 
is convenient for comparison with the analysis of 
Ref. 11. 


Theory 

In this section, a brief review will be pre- 
sented of a "steady" state finite element theory and 
a transient finite difference theory. In the 
following section of this report, these theories 
will be compared to the experimental data. 

Finite Element Tneory 

For the case of a sheared mean flow in a rec- 
tjangular duct, Astley and Eversman^^ nave devel- 
oped a finite element method based on weighted 
residuals using an eight node isoparametric element, 
snown in Fig. 4. 

The acoustic equations are described by the 
linearized perturoed energy and momentum equations 
for isentropic flow. The perturbed quantities are 
assumed to vary with time as e''“**^*, and the 
( ‘eory thereby is called "steady" state theory, 
Velocities and pressures are nondimensionalized 
using reference 'cales, cS and ppcj'^, 
respectively. 1 ■» following equations result: 

31) 3U 

(x-mom.) ixu + H + u ^ + V ^ 


1 Dp 1 ®*^0 


(2) 


(y-mom.) iku + + V^ + 




+ L.£_ 

Pq y 


1 ^Pq 
YPoPp 3y 
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(3) 


Dp Dp 

(energy) Ikp + U + v + u -r” v 

+ Jp t p + p!Sa t P S’] . 0 (4) 

’ ro Dx ♦'o Dy *^Dx ^ Dy J '' ' ' 

The complete details for transforming these equa- 
tions Into the finite element formulation can be 
found In Ref. 11. 

The speaker, see Fig. 1, sends acoustic pres- 
sure waves Into the duct which are either reflected 
or transmitted by the variable area duct restric- 
tion. Accurate source modelling and representation 
of the duct termination requires special considera- 
tion, Acoustic mode reflection at the Inlet to the 
duct nonuniformity and the transmission without 
reflection at the outlet of the nonuniformity Is 
modelled by matching the finite element solution In 
the Interior of the nonuniformity to analytical or 
finite element eigen function expansions in the uni- 
form inlet and outlet ducts. This permits a multi- 
modal representation accounting for reflection and 
mode conversion by the nonuniformUy. This approach 
was first introduced in the acoustic problem by 

Eversmaq, ,Cook and Beckemeyer^'* In connection with 
a modified GalerKln solution for transmission in 
nonuniform ducts without flow and was subsequently 
used in a similar context in the presence of 

fl owl 1.1b. 


The boundary condition at the upper and lower 
hard walls requires that the normal acoustic veloc- 
ity be zero. 

Transient Finite Difference Theory 

Will tei^ has numerically mapped the variable 
area duct into the rectangular geometry shown in 
Fig. b. The wave equation In the physical xy plane 

was transformed to the ?,n plane 


(b) 


Where J,o,a,7]_,o and t are functions of x,y,z 
and n and are defined in Ref. 12. The mapping 
trom the physical to the transformed C.n plane was 


performed numerically using a technique developed 
by Thompson (Refs. Id and 1?), Equation (6) was 
tnen rewritten in finite difference form and solved. 


The entrance and exit conditions used in the 
transient finite difference theory were 


P(o.y.t) , 


(7) 


and an exU iiDpeaance value of 

2%/L.1-pocS (8) 
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itoth Eqs. (7) and (8) represent approximations to the 
entrani‘4 and exit conditions employed by Astley and 
Eversman (Ref. 11) discussed earlier. At the entrance 
and exit of the variable area section, higher order 
acoustic inodes v^lll generally yield a non-uniform 
pressure profile and a noivuniform exit Impedance. 

The degree which Eqs. (7) and (8) approximate the true 
entrance conditions can be Improved by adding both 
straight entrance and exit sections. In these sec- 
tions. the hlaher order modes will decay such that 
Eqs. (7) and (8) will more closely approximate the 
true boundary conditions. 

Olscusslon of Results 


In this section, experimental results are dis- 
cussed In the following categories; rms axial pres- 
sure profiles, transverse pressure profiles and phase 
angle. 

RMS Axial Pressure Profiles 


The normalized rms pressures from Table I are 
plotted against axial distance In Fig. 6 and compared 
to the theories of Refs. 11 and 12. As seen in 
Fig. 6, the Astley-Eversman theory Is In excellent 
agreement with the data. In their analysis, the 
Incident pressure wave had a value of 1; thus, the 1.4 
value at x « 0 Indicates a reflected wave of 0.4 in 
magnitude. 

The White transient finite difference theory 
snown by the dashed line In Fig. 6 is also in reason- 
ably good agreement with the data. Some deviation of 
the pressure at x*/l* equals 1.0 Is seen. 

Transverse Pressure Profile 


Figure 7 shows a comparison between the measured 
and calculated transverse pressure profiles. Both 
theories predict very small differences in pressure 
oetween the flat (y = o) and the curved surface 
(y = 1). The scatter In the data appear large because 
the profile is nearly flat; thus the pressure differ- 
ences In Fig. 6 are very small. Also, the pressure 
at X P 0.3 on the flat surface had to be estimated 
from Fig. 6 because a microphone was not at this exact 
location on the flat surface. Nevertheless, both 
theories predict tne proper trends of the data and are 
in reasonable agreement in magnitude with the test 
data. 

Phase Angle 

Figure 8 shows a comparison between the measured 
and calculated value of the phase angle as a function 
of axial distance. Again both theories are in reason- 
able agreement with tne data. In this case, the 
transient finite difference theory gives slightly 
better agreement at x*/L* of 1. Overall, consider- 
ing Figs. 7 and 8, both theories appear to perform 
equally well in correlating the data. 

Standing Wave Pattern 

The duct construction In the present experiment 
produced a reflection which induced a significant 
standing wave pattern ahead of the test section. 

Also, any reflection from the duct termination will 
generate standing waves downstream of the test sec- 
tion. To determin'' this wave pattern, rms pressure 
measurements were taken along the flat wall both up- 
stream and downstream of the test section. These mea- 
surements are recorded In Table III and plotted In 
rig. y. As seen in Fig. y, tne finite element theory 


and the transient finite difference theory are In 
excellent agreement with the data upstream of the test 
section. Since the theory assumes that there Is no 
reflected energy at the exit, a downstream standing 
wave pattern cannot, of course, arise. 

The standing wave ratio S downstream of the 
test section was 


S 


0.905 

or 


1.09 


(9) 


Therefore, assuming the Incident and reflected waves 
react In normal Incidence at the duct exit, the sound 
power reflection coefficient Is (Ref, 18, p. 301) 

*r " “ 0.0018b (10) 

This slight amount of reflected sound from the duct 
exit probably accounts for part of the deviation be- 
tween experiment and theory which assumes a non- 
reflecting exit termination. 

With the exception of the Astley-Eversman modal 
coupling finite element formulation, all other numeri- 
cal analyses specify an entrance pressure or velocity 
and an approximate exit Impedance, white's analysis, 
which uses Eqs. (7) and (8) for the entrance and exit 
conditions, is typical of these theories. Direct pre- 
diction of the transmitted and reflected energy Is not 
generally possible In these theories. For low fre- 
quency plane wave propagation (Ref. 18, p. 300, 

Ref. 19), however, an estimate of the reflection 
coefficient can be found by adding a straight section 
of duct ahead of the area variations and measuring the 
standing wave pattern In the region where higher order 
modes have decayed. Higher order modes will be gen- 
erated at the interface between the straight and vari- 
able area portion of the duct. 

In contrast, the Astley-Eversman finite element 
formulation predicts the transmitted and reflected 
energy by analytical Including a straight entrance and 
exit section; therefore, their analysis does not re- 
quire elements In these sections. Application of the 
Astley-Eversman theory to only the contracting portion 
of the duct yielded a reflection coefficient of 
0.264. This coefficient can now oe compared to the 
reflection coefficient obtained from the measured and 
predicted standing wave pattern shown in Fig. 9. 

From Fig. 9, the standing wave ratio Is 


<; 1.52 

^ “ 0.49 


3.120 


( 11 ) 


Since only plane waves propagate in the straight por- 
tion of the duct at w 0.172, the power reflection 
coefficient can be calculated directly (Ref. 18, 

£q. 8.49). 


„ - ( 3.102 - 1 )^^ - 

(3.102 + 1)*^ 


0.263 


( 12 ) 


Which is in excellent agreement with the value pre- 
dicted by the Astley-Eversman modal coupling theory. 
Thei^efore, in a general situation where it is desired 
to estimate the reflected or transmitted acoustic 
energy, the eigenvalue formulation of Astley and 
Eversman for the duct entrance and exit condition 
would be the appropriate algorithm to use. There 
could be some difficulty, however, in estimating the 
modal content of the incident wave if higher order 
modes can propagate. 
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Conclusion 

Experimental data have been presented for sound 
propagation in a simulated infinite hard wall duct 
with a large change in duct cross sectional area, The 
data is conveniently tabulated for further use. The 
"steady** state finite element theory of Astley and 
Eversman and the transient finite difference theory of 
White are in good agreement with the data for both the 
axial and transverse pressure profiles and the axial 
phase angle. Therefore, numerical finite difference 
and finite element theories appear to be ideally 
suited for handling duct propagation problems which 
encounter large axial gradients in acoustic parameters. 

The measured energy reflection coefficient agrees 
with the values from the Astley-Eversman modal coupl- 
ing model. Therefore, in a general multi-modal situa- 
tion where it is desired to estimate the reflected or 
transmitted acoustic energy, the eigenvalue formula- 
tion of Astley and Eversman for the duct entrance and 
exit conditions would be the approximate algorithm to 
use. Definition of the mode content of the incident 
wave might represent a problem, however. 
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TABLE I. - AXIAL PRESSURE MEASUREMENTS 
Straight Mall 


Data Set I. Prms j 

Data Set 

• i 1 Prits 

1 

2 

3 

Average 

1 

2 

3 

Average 

1.395 

1.390 

1.418 

1.40 

1.405 

1.403" 

1.395 

1.40 

1.556 

1.473 

1.508 

1.513 

1.443 

1.468 

1.493 

1.468 

1.490 

1.508 

1.488 

1.495 

1.500 

1.480 

1.500 

1.493 

1.397 

1.385 

1.387 

1.390 

1.405 

1.380 

1.400 

1.395 

1.252 

1.234 

1.239 

1.242 

1.214 

1,217 

1.237 

1.222 

1.066 

1.069 

1.046 

1.061 

1.049 

1.066 

1.061 

1.059 

.918 

.918 

.916 

.918 

.941 

.933 

.941 

,938 

.835 

.841 

.868 

.848 

.851 

.851 

.863 

.856 

.810 

.825 

.825 

.820 

.805 

.815 

.823 

.815 


Curved Wall 


MIC i 
position 
X*/L* 


Odtd S6t If 

Data Set 

II, Prms 

1 

2 

3 

Average 

1 

2 

3 

Average 

-0.126 

1.405 

1.443 

1.440 

1.43 

1.387 

1.382 

1.423 

1.397 

.3 

1.573 

1.535 

1.578 

1.536 

1.490 

1.553 

1.551 

1.530 

.5 

1.259 

1.254 

1.227 

1.247 

1.219 

1.239 

1.295 

1.219 

.7 

.956 

„938 

1.006 

.966 

.936 

.911 

.951 

.933 

1.125 

.876 

.883 

.886 

.881 

.848 

.863 

.891 

• 368 


TABLE II. - PHASE ANGLE DATA 
Straight wall 


MIC t 

Data Set I, ' 

Data Set II, P 

position 

X*/L* 

Bi 

2 

3 

Average 

B 

2 

3 

Average 

0.0 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

.125 

5.70 

6.5 

6.2 

6.13 

6.1 

6.3 

6.2 

6,2 

.250 

12.6 

13.0 

12.9 

12.83 

13.0 

13.0 

13.0 

13.0 

.375 

21.2 

22.0 

21.3 

21.5 

21.4 

21.4 

21.4 

21.4 

.50 

32.4 

33.3 

32.8 

32.83 

32.8 

32.8 

32.9 

32.83 

.626 

47.9 

48.6 

48.4 ! 

48.3 

48.1 

48.1 

48.3 

48.17 

.750 

66.0 

66,6 

66.6 

66.4 

66.6 

66.7 

66.7 

66.67 

.875 

86.4 

86.9 

86.1 

86.47 

86.4 

86.5 

86.5 

86.47 

1.0 

105.6 

106.6 

106.8 , 

106.33 j 

106.4 

106.4 

106.4 

106.4 


Curved Wall 


MIC t 


Data Set I, p 



Data Set II, ^ 

position 

X*/L* 

1 

2 

3 

Average 

1 

2 

3 

Average 

-0.125 

-4 ,47 

-4.27 

-4.87 

,4.64 

-5.10 

mm 

-5.D 

-5.13 

.3 

10,43 

10.93 

10.43 

10.6 

10.1 

mimm, 

10.1 

10.27 

.5 

32.83 

32.63 

32.43 

32.63 

32.3 j 

32.5 

32.3 

32.37 

.7 

71.23 

72.13 

72.73 

72.03 

71.7 ' 

71.4 

71.5 

71.53 

1.125 

119.73 

119.83 

120.13 

119.9 

119.2 ! 

119.5 

119.4 

119.37 
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Figure 1. ~ No flow acoustic duct test section and instrumentation. 
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Figure 2 . • Variable area test piece with 3 microphone locations 
(microphone plugs not shown). 


Figure 3. - Microphone test holder plate mounted above the variable area 
test piece. 
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Figure 4. - Finite element discretization of variable ares durt using theory 
of reference 11. 




(bl MAPPED PLANE. 

Figure 5. - Mapping for a variable area duct 












